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ABSTRACT 

We compile a new sample of 1 15 nearby, non-Seyfert galaxies spanning a wide range of star 
' formation activities, from starburst to nearly dormant, based on ultraviolet observations with 

various satellites. We combine these observations with infrared observations to study the re- 
lation between ratio of total far-infrared to ultraviolet luminosities and ultraviolet spectral 
slope (the 'IRX-UV relation). We show that, at fixed ultraviolet spectral slope, quiescent 
star-forming galaxies in our sample have systematically lower ratio of total far-infrared to 
ultraviolet luminosities than starburst galaxies. The strengths of spectral indices sensitive to 
| star formation history, such as the 4000 A spectral discontinuity and the Ha emission equiv- 

C^) • alent width, correlate well with distance from the mean relation for starburst galaxies in the 

IRX-UV diagram, while there is little or no correlation between the dust-sensitive Ha/H/3 
ratio and this distance. This is strong observational evidence that the star formation history 
is relevant to the 'second parameter' affecting the IRX-UV relati on. We show that th ese re- 
sults can be understood in the framework of the simple model of Chariot & Fall 2000 for the 
transfer of starlight through the interstellar medium in galaxies. We confirm that, for starburst 
galaxies, the tight IRX-UV relation can be understood most simply as a sequence in overall 
dust content. In addition, we find that the broadening of the relation for quiescent star-forming 
galaxies can be understood most simply as a sequence in the ratio of present to past-averaged 
star formation rates. We use a library of Monte Carlo realizations of galaxies with different 
r> | star formation histories and dust contents to quantify the accuracy to which the ultraviolet 

^ ■ attenuation Apuv of a galaxy can be estimated from either the ratio of far-infrared to ultravi- 

olet luminosities or the ultraviolet spectral slope. We provide simple formulae for estimating 
Apuy as a function of either of these observational quantities and show that the accuracy of 
these estimates can be significantly improved if some constraints are available on the ratio of 
present to past-averaged star formation rates. 
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1 INTRODUCTION 

Much of what we know about the universe at high redshift arises 
from observations of the rest-frame ultraviolet emission of distant 
galaxies that is accessible from the ground at optical and near- 
infrared wavelengths. For this reason, significant effort has been 
devoted over the past several years to improve our understand- 
ing of the star formation history and dust attenuation in galax- 
ies drawn from ultraviolet surveys. Recent progress in this area 
has come mainly from studies of a sample of 57 nearby starburst 
galaxies, for which a wide range of observations were compiled by 
M eurer, Heckman & Calzetti il999l) . In particular, for these galax- 



ies there is a remarkably tight correlation between ratio of total 
far-infrared to ultraviolet luminosities, -Ldust/I/Fuv, and ultravi- 
olet spectral slope, (3, which is referred to as the 'infrared excess 
(IRX)-ultraviolet (UV)' relation (here (3 is defined by a power- 
law fit of the form fx oc A^ to the spe c trum at wave lengths 
1300 < A < 2600 A: ICalzetti et aljfl994h . IMeurer et"al] <1999h 
find that the IRX-UV relation allows reliable estimates of the at- 
tenuation by dust at ultraviolet wavelengths, Afuv, from either 
^dust/^FUv or j3, This relation, if universal, represents our best 
hope to estimate the attenuation by dust and hence the star forma- 
tion rates of high-redshift galaxies, for which the rest-frame optical 
and infrared emission cannot be observed at present. 
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The main uncertainty affecting the IRX-UV relation is that 
it has been established only for starburst galaxies. Recently, iBelll 
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(2002) has compiled a small sample of nearby, more quiescent 
galaxies drawn from ultraviolet observations with various satellites 
(here we use the term 'quiescent' to describe the broad class of 
star-forming galaxies between active starburst and dormant early- 
type galaxies). His results suggest that quiescent galaxies deviate 
from the IRX-UV relation for starburst galaxies, in the sense that 
quiescent galaxies tend to have 'redder' ultra violet spectra a t fixed 
idust / '£fuv ■ This implies that the simple IMeurer et all 1 19991) 
recipe to estimate Afuv from either L d us t /£fuv or /? may not 
hold for all galaxy types. The results of Bell 1 2002) raise several 
questions. Can star formation history alone account for the distri- 
bution of galaxies in the I RX-UV diagram? Or d o the galaxies have 
different dust properties? Charl ot& Falll l2000l) have developed a 
simple dust model, based on an idealized prescription of the main 
features of the (clumpy) interstellar medium (ISM), which allows 
one to interpret the IRX-UV relation and the emission-line prop- 
erties of starburst galaxies in terms of star formation rate and dust 
content. This model accounts for the different attenuation affect- 
ing young and old stars in a galaxy, because of the dispersal of 
the (g iant molecular) clouds in which stars form. Using this model, 
Chari ot & FaTilfeOOOl) show that the tight IRX-UV relation for star- 
burst galaxies can be understood most simply as a sequence in the 
overall dust content. Clearly, to gain insight into the origin of the 
IRX-UV relation for different types of galaxies, we require further 
observations and modelling of not only the ultraviolet and infrared 
luminosities, but also other observable quantities with distinct de- 
pendences on dust content and star formation history. 

In this paper, we compile a new sample of 115 nearby, non- 
Seyfert galaxies spanning a wide range of star formation activities, 
from starburst to nearly dormant, to explore the physical parameters 
driving the location of a galaxy in the IRX-UV diagram. Our sam- 
ple is drawn from ultraviolet observations with various satellites. 
We combine these observations with infrared observations and, for 
a subset of the galaxies, with optical spectroscopy from various 
sources. We examine, in particular, the 4000 A spectral disconti- 
nuity and the Ha emission equivalent width, which are good indi- 
cators of star formation activity, and the Ha/H/3 ratio, which is a 
good indicator of attenuation by dust. Based on these data alone, it 
appears that the broadening of the IRX-UV relation for quiescent 
star-forming galaxies is most likely driven by star formation his- 
tory. We combine the simple dust model of Chariot & Fall 1 2000 ) 
with the population synthesis code of Bruzual & Chariot 1 2003) 
and show that this model can account for the ultraviolet and far- 
infrared properties of the galaxies in our sample. The model al- 
lows us to understand why the ratio of present to past-averaged star 
formation rates is the most likely 'second' parameter that, in com- 
bination with dust content, determines the position of a galaxy in 
the IRX-UV diagram. In Section 2 below, we compile our galaxy 
sample and investigate the origin of the IRX-UV relation from an 
observational viewpoint. In Section 3, we interpret the IRX-UV re- 
lation using the model. We explore the implications of our results 
for estimates of the ultraviolet attenuation in galaxies in Section 4. 
Our conclusions are summarized in Section 5. 



2 THE ORIGIN OF THE IRX-UV RELATION: 
OBSERVATIONAL INSIGHT 

In this section, we re-examine the IRX-UV relation from an ob- 
servational viewpoint. We combine ultraviolet observations from 
various satellites with observations from the Infrared Astronomical 
Satellite (IRAS) to compile a new sample of nearby, non-Seyfert 



galaxies spanning a wide range of star formation activities in the 
IRX-UV diagram. We explore how the position of a galaxy in 
this diagram depends on observable properties with distinct depen- 
dences on dust content and star formation history. 

2.1 The data 

We compile a new sample of nearby, non-Seyfert galaxies span- 
ning a wide range of star formation activities, using available ob- 
servations with the International Ultraviolet Explorer (IUE), the 
Orbiting Astronomical Observatory (OAO) and the Astronomical 
Netherlands Satellite (ANS). We first extract all (non-Seyfert) spi- 
ral, elliptical, irregular and emission-line galaxies (classes 80, 81, 
82 and 88) from the IUE Newly-Extracted Spectra (INES version 
3.0) archive. The spectra in this database extend over two wave- 
length domains, corresponding to the far- (1 150-1980 A) and near- 
ultraviolet (1850-3350 A) spectrographs on-board IUE. We re- 
ject galaxies with only near-ultraviolet spectra, for which the ul- 
traviolet spectral slope cannot be measured reliably. When both 
the far- and near-ultraviolet spectra are available for a galaxy, we 
merge them to form a single ultraviolet spectrum encompassing the 
range from 1150 to 3350 A (the spectra usually agree to within 
a few percent in the overlapping region around 1900-1950 A). 
We remove foreground Ga lactic extin ction using the dust maps 
of Schleael, Finkbeiner & Davis 1 1998) and the Galactic extinction 
curve of Cardelli, Clavton & Mathis 1 1989). Furthermore, we co- 
add the spectra of galaxies observed multiple times to increase the 
signal-to-noise ratio. For the ultraviolet spectral slope to be mea- 
sured with reasonable accuracy, a minimum signal-to-noise ratio 
per pixel of 3 is required (see Appendix). A total of 110 galaxies 
with both far- and near-ultraviolet spectra and 185 galaxies with 
only far-ultraviolet spectra satisfy this cut in signal-to-noise ratio. 

We require two additional conditions to include galaxies in our 
final sample: that optical photometry be available to estimate the ul- 
traviolet aperture correction, and that IRAS 25, 60 and 100 /im ob- 
servations be available to estimate the total luminosity reradiated 
by dust. The above sample of 295 galaxies includes the 57 star- 
bu rst galaxies with opt ical angular diameters less than 4' studied 
by Meurer et al. 1 1999), for which most of the concentrated ultra- 
violet emission could be observed within the 20" x 10" IUE aper- 
ture. IRAS 25, 60 and 100 fim observations are available for 50 of 
the Meurer et al. galaxies, which we include in our final sample. 
For other galaxies, reasonable ultraviolet aperture corrections may 
be derived as a function of the effective o ptical radius and mor- 
phological type iRifatto. Lonao & Capaccioli 1995b). Morpholog- 
ical types are known for IUE galaxies, but effective optical radii 
must be obtained from independent observations. We thus cross- 
correlate th e above sam ple of 238 IUE galaxies not included in 
the M eurer et alj 1 19991) sample with th e Data Release One of the 
Sloan Digital Sky Survey (SDSS DR1 : lAbazaiian et aljEooll) . In 
this way, we obtain effective g-band radii for 41 galaxies (with 
typical, i.e. median r e fi ~ 9"), for which we compute ultraviolet 
aperture corrections using the prescription o flRifatto et aljfl995bl) . 
The corrections vary from galaxy to galaxy, with a mean logarith- 
mic ratio of corrected-to-observed ultraviolet fluxes of 0.2 ± 0.2. 
Only 33 of these galaxies have available IRAS 25, 60 and 100 fim 
observations and can be included in our final sample. These typi- 
cally have morphological types earlier than the starburst galaxies 
of the IMeurer et alj 1 1999) (Sa-Sc versus Sc-Im, although some 
starburst galaxies are also present) and similar absolute B-band 
magnitudes, —16 < Mb < —22 (we adopt a Hubble constant 
H = TOkms^Mpc" -1 ). 
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We supplement this sample with (non-Seyfert) galaxies for 
which multi-band ultraviolet imaging with ANS and OAO is avail - 
able from the compilation o flRifatto. Longo & Canacciolil ll995al) . 
In the spectral range of interest to us, the galaxies observed with 
ANS were imaged at 1550, 1800, 2200 and 2500 A through a fixed 
aperture of 2.5' x 2.5'. Those observed with OAO were imaged at 
1550, 1910 and 2460 A through a fixed circular apertu re of 10' di- 
ameter. Of the 95 galaxies observed with ANS in the Rifatt o et alJ 
1 1995a) compilation, 34 have available IRAS 25, 60 and 100 pm ob- 
servations. We reject 5 galaxies with optical diameters larger than 
9', for which the ultraviolet emission does not appear to be clearly 
concentrated within the ANS aperture, and 2 galaxies with larger- 
aperture imaging available from OAO (see below). The remaining 
27 galaxies have typical morphological types Sb-Sd and absolute 
B-band magnitudes similar to those of the galaxies in the IUE sam- 
ple, -17 < M B < -22. Of the 36 galaxies with OAO imaging, 11 
have available IRAS 25, 60 and 100 /xm observations. We reject 6 
galaxies with optical diameters larger than 12'. We thus retain only 
5 galaxies with OAO imaging in our final sample, two Sc galaxies 
with Mb ~ —16, two Sd galaxies with Mb ~ —20 and an Im 
galaxy with Mb ~ —19. Our final sample contains 115 galaxies 
with both ultraviolet and infrared data, of which 83 were observed 
with IUE, 27 with ANS and 5 with OAO. About ha lf are starburst 
galaxies (most belonging to the lMeurer et all 19991 sample), while 
the other half are earlier-type, more quiescent galaxies spanning a 
similar range of absolute B-band magnitudes. 

We comp ute ultraviolet spect ral slopes for all the galaxies 
in our sample. ICalzetti et alJ 1 19941) define the ultraviolet spectral 
slope (5 by a power-law fit of the form fx oc A* 3 (where f\ is the 
flux per unit wavelength) to the spectrum in ten continuum bands 
in the range 1268 ^ A ^ 2580 A. We cannot use this definition 
here, because it does not apply to galaxies for which only multi- 
band ultraviolet imaging is available. Instead, we adopt the follow- 
ing definition of the ultraviolet spectral slope for all the galaxies in 
our sample, 



/3 G LX = lQg /FUV ~ lo S /NUV 

log Apuv — log Anuv ' 



(1) 



where Afuv = 1520 A and Anuv = 2310 A are the effective 
wavelengths of the far- (1350-1800 A) and near-ultraviolet (1800- 
3000 A) filters on board the Galaxy Evolution Explorer (GALEX, 



iMartin et all2003l) . and / PUV and / NUV are the mean flux densi- 
ties (per unit wavelength) through these filters. In the Appendix, we 
use a reference sample of 1 10 IUE spectra to calibrate the conver- 
sions between /3glx and the ultraviolet spectral slopes estimated 
from observations with different satellites. We show that, in partic- 
ular, ultraviolet spectral slopes estimated from ANS and OAO imag- 
ing observations can be converted into /3glx with an rms uncer- 
tainty of only a ~ 0.2. This is similar to the uncertainty introduced 
in /3glx by the lack of near-ultraviolet spectra for some IUE galax- 
ies in our sample. For all the galaxies in our sample, we compute 
the ultraviolet flux as Ffuv = <^fuv/ftjv> as described in the 
Appendix. 

To estimate the total far-infrared flux Fd ust from the observed 
IRAS flux densities /i/(25 M-m), ^ (60 urn) and /WlOO iim), we 
appeal to the recent prescription of Dale & Helou (2002). This is 
based on model fits to the spectra of a large sample of nearby, nor- 
mal (i.e. non-Seyfert) star-forming galaxies across the whole wave- 
length range from 3 to 850 /im. We use equation (5) of Dale & 
Helou to compute Fdust CFtir in their notation) from f v (25 pm), 
/„(60 pm) and /„(100 \im). For reference, for the IUE starburst 
galaxies in our sample, the values of Fdust computed in this way are 
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Figure 1. Ratio of far-infrared to ultraviolet luminosities plotted against 
ultraviolet spectral slope (the 'IRX-UV diagram'). The data points are from 
the sample compiled in Section 2.1 (median measurement errors for each 
subsample are indicated in the lower right corner). The solid line shows the 
mean IRX-UV relation for the 50 starburst galaxies in this sample (equa- 
tion0. 



typically 50% larger than those estimated by M eurer et alJ 1 19991) 
from only /„(60 and f^flOO \im) usi ng the older prescrip- 
tion of Helou et al. (1988; see also lCalzetti et al]l2000h . Neglect- 
ing possible anisotropies (Section 4), we equate the ratio of total 
far-infrared to ultraviolet fluxes Fd us t / i*Fuv to the corresponding 
luminosity ratio I/dust /£fuv- 



2.2 The observational IRX-UV diagram 

Fig.Qshows Ldust /Lfuv as a function of /3glx for the 115 galax- 
ies in our sample. Different symbols distinguish the 50 starburst 
galaxies studied by Meurer et al. (1999, solid dots) and the (typ- 
ically more quiescent) galaxies observed with IUE+SDSS (stars), 
ANS (open triangles) and OAO (filled triangles). To a first, rough 
approximation, £d us t / Lfuv is a function of /3glx for the whole 
sample. As found previously by Meurer et al. (1995, 1999), the two 
quantities are tightly correlated for starburst galaxies, for which the 
typical scatter at fixed /Sglx amounts to a factor of less than 2 
in iduat/^FUv- The main novelty introduced by the inclusion of 
more quiescent galaxies in our sample is an increase of the scat- 
ter to a factor of more than 3 toward lower Ldust / L-fuv values at 
fixed /3glx- The effect is real, as the observational scatter is much 
larger than the typical (i.e. median) measurement errors (which are 
shown in the lower right corner of Fig.0. We note that the shape of 
the IRX-UV relation does not appear to be affected by the inclusion 
of quiescent galaxies in our sample. 

What are the main physical parameters driving the relation be- 
tween ultraviolet spectral slope and ratio of far-infrared to ultravi- 
olet luminosities for star-forming galaxies? ICharlot & Falll 1 2000) 
show that, for starburst galaxies, the tight relation between /3glx 

1 It is worth noting that no galaxy in Fig. lies substantially above the 
IRX-UV relation for starburst galaxies. This is not unexpected, since our 
sample of galaxies drawn from ultraviolet surveys (with typical infrared 
luminosities Ldust ~ 10 10 Lq) does not include 'luminous and ultralu- 
minous infrared galaxies', which tend to have larger idust/^FUV than 
starburst galaxies at fixed /3glx (presumably because of the heavy obscu- 
ration of the sources of ultraviolet radiation; Goldader et al. 2002). 
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and Ldust /£fuv can be understood most simply as a sequence in 
the overall dust content of the galaxies. The inclusion of quiesent 
star-forming galaxies in Fig.Q reveals that a 'second parameter' 
tends to broaden the relation. We may begin by identifying this 
second parameter as the perpendicular distance fr om the mean 
IRX-U V relation for starburst galaxies. Following Meur er et alj 
( 1999), we define the mean relation for starburst galaxies by a least- 
squares fit of the type 

L dust /L FUV = 10( 2 - 10+0 - 85 ^ GLX ) - 0.95 (2) 

to the observations of the 50 starburst galaxies in our sample. The 
result is shown by the solid line in Fig.Q We can compute the 
perpendicular (i.e. shortest) distance dp from each point in the 
IRX-UV diagram to this line. We adopt the convention of nega- 
tive dp for galaxies with Ldust /Lfuv lower than the mean value 
for starburst galaxies at fixed /3glx (see Fig.0. 

We expect the second parameter affecting the IRX-UV 
relation to be connected to star formation history, since the 
relation becomes broader when quiescent star-forming galaxies 
are considered in addition to starburst galaxies. It is thus natural 
to explore how observed quantities with different dependences on 
the star formation history correlate with dp. We first consider the 
4000 A spectral discontinuity index D n (4000), defined as the ratio 
of the average flux densiti es in the narro w bands 3850-3950 A 
and 4000-4100 A {Bal oghetalJ Il999l) . This index depends 
somewhat on metallicity but correlates more with the ratio of 
present to past-averaged star formation rates in galaxies, and 
hence, it is a valuable i ndicator of the history of star formation 
iBrinchmann et al]|2003l) . Also, the narrowness of the bands used 
to compute D n (4000) makes it largely insensitive to attenuation by 
dust. Medium-resolution spectra (3-9 A FWHM) of 64 galaxies in 
our sampl e (35 starburst an d 29 quiescent s tar-formin g) are avail - 
able from|Kennicuttl ll992T).lMcOuade. Calzetti & Kinnevl Il995l>. 
IStorchi-Bergrrra^rrTKm^ Ji995Tk iKong&CnensI 
120021) and lAbazaiian et alJ <2003l) . Fig.^ shows D n (4000) 
measured from these spectra as a function of the perpendicular 
distance dp from the mean starburst IRX-UV relation. The 
Spearman rank correlation coefficient is r s = —0.61, indicating 
that the correlation between D n (4000) and dp is significant at the 
5 a level for this sample size. This result is remarkable, as it is 
direct observational evidence that star formation history is relevant 
to the second parameter affecting the IRX-UV relation. 

Another observable quantity that is sensitive to the ra- 
tio of present to past-averaged star forma tion rates is the 
Ha emission equivalent width E W(Hq) iKennicuttl Il983t 
iKennicutt. Tamblvn & Cong don 1994). This is the ratio of the Ha 
emission-line luminosity produced by young massive stars to the 
continuum luminosity per unit wavelength produced by older stars 
at Ahq = 6563 A. In Fig.[2|}, we show EW(Ha) as a function of 
dp for 59 galaxies (31 starburst and 28 quiescent star-forming) in 
our sample (we adopt the convention of positive equivalent width 
for emission). The Ha emission equivalent widths were measured 
from the same optical spectra as above following the procedure out- 
lined in section 2 of Brinchmann et al. (2003), which allows accu- 
rate correction for stellar absorption. The Spearman rank correla- 
tion coefficient in this case is r s — 0.48, indicating that the cor- 
relation between EW(Ha) and dp is significant at the 4 a level for 
the sample size. We note that dust can reduce EW(Hq) through the 
absorption of ionizing photons. However, the effect is modest for 
quiescent star-forming galaxies, typically only 2 0-30 per cent, an d 
cannot account for the trend seen in Fig.|2l) ^harlot et alll2002h . 
Taken together, Figs|2^ and b show that two different tracers of the 



ratio of present to past-averaged star formation rates, D„ (4000) 
and EW(Hq), correlate with dp. This is strong evidence that the 
ratio of present to past-averaged star formation rates could be the 
second parameter affecting the IRX-UV relation. 

In an attempt to further understand the relative influence of 
star formation history and attenuation by dust on the IRX-UV 
relation, we plot other properties of the galaxies against dp in 
Figs|2j; and d: the ratio of the IRAS /„(60 ]xm) and /„(100 u.m) 
flux densities and the ratio of the Ha and Hf3 emission line lumi- 
nosities. The /„(60/im)//„(100/im) ratio, available for all 115 
galaxies in our sample, depends on the relative contributions by 
young and old stars to dust heating in a galaxy (e.g.. lHeloull986l) . 
Hence it is an indirect tracer of the star formation history that 
also depends on the dust content (and potentially the spatial dis- 
tribution of dust relative to stars). The Ha/H/3 ratio, available for 
48 galaxies (29 starburst and 19 quiescent star-forming) from the 
same optical spectra as above, is strongly sensitive to attenuation 
by dust but not to the past history of star formation. The Spear- 
man rank correlation coefficients are r s = +0.43 for the rela- 
tion between /„(60 /jm)//„ (100 /im) and dp and r s = 0.12 
for that between Lna/Lnp and dp. For the different numbers of 
galaxies in each case, this implies that the correlation between 
/„(60/im)//„(100/im) and dp is significant at the 4 a level, 
while there appears to be little or no correlation between Lh[*/Lh/3 
and dp. These findings reinforce our contention that the second pa- 
rameter affecting the IRX-UV relation is primarily star formation 
history rather than dust content. 

Based on the observational evidence above, we conclude that, 
while dust content is the primary driver of the IRX-UV relation 
for starburst galaxies, the broadening of the relation for quiescent 
star-forming galaxies is most likely driven by star formation his- 
tory, and in particular the ratio of present to past-averaged star for- 
mation rates. In the next section, we show that this observational 
result can be understood in the framework of a simple model for 
the production of starlight and its transfer through the interstellar 
medium (ISM) of galaxies. 



3 THE ORIGIN OF THE IRX-UV RELATION: 
MODELLING 

In this section, we explore the physical origin of the IRX-UV re- 
lation from a more theoretical standpoint. We appeal to a simple 
model for the production of starlight and its transfer through the in- 
terstellar mediu m in galaxies. Th e model is based on a combination 
of thelBruzual & Charlotl {2003) population synthesis code and the 
ICharlot & FallH2000l) prescription for the absorption of starlight by 
dust. The Chariot & Fall model relies on an idealized description 
of the main features of the ISM. Stars are assumed to form in in- 
terstellar 'birth clouds' (i.e. giant molecular clouds). After 10 7 yr, 
young stars disrupt their birth clouds and migrate into the 'ambient 
ISM'. The 'effective absorption' curve describing the attenuation 
of photons emitted in all directions by stars of age t' in a galaxy is 
given by the simple formula 

/ tv (A/5500 A) -0 ' 7 , fort's; 10 7 yr, 
T\(t ) = < } (_ 7 (3) 

[ [it v (A/5500 A) , for t' > 10 7 yr, 

where fv is the total effective V-band optical depth seen by young 
stars. The adjustable parameter n defines the fraction of the total 
effective absorption optical depth contributed by the ambient in- 
terstellar medium (n ~ 1/3 on average, with substantial scatter). 
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Figure 2. (a) spectral discontinuity at 4000 A, (b) Ha emission equivalent width, (c) ratio of IRAS /„(60 |jm) and /„(100 |am) flux densities and (d) Ha/H/3 
ratio plotted as a function of perpendicular distance from the mean starburst IRX-UV relation, for the subsets of galaxies for which these observed quantities 
are available in the sample of Fig.0 The different symbols have the same meaning as in Fig.fTl 



In addition, a fraction / m 0.1 of the effective absorption optical 
depth (1 — f£}Tv of the birth clouds is assumed to arise from the 
Hll regions ionized by young stars in the inner parts of these clouds. 
We note that, in the case of optically thick birth clouds, the ultra- 
violet radiation escaping from a galaxy is produced by stars older 
than ~ 10 7 yr, while the total far-infrared luminosity is produced 
by s tars of all ages. 

ICharlot & Falll l2000h used their model to interpret the rela- 
tions between ultraviolet spectral slope, ratio of far-infrared to ul- 
traviolet luminosities and Ha/H/3 ratio for the starburst galaxies 
in the Meu rer et alJ (1999) sample (Section 2). They showed that 
the finite lifetime of the birth clouds is a key ingredient for re- 
solving the apparent discrepancy between the attenuation of line 
and continuum radiation in these galaxies. They also showed that 
the IRX-UV relation can be understood most simply as a se- 
quenc e in the overall dust content of the galaxies. For simplicity, 
Charl ot"& Falll j2000) limited their study to models with constant 
star formation rates and young ages, as appropriate to interpret ob- 
servations of starburst galaxies. The results of Section 2 above sug- 
gest that star formation history plays an important role in defining 
the IRX-UV relation for more quiescent galaxies. To investigate the 
dependence of the IRX-UV relation on star formation history, we 
assume that the star formation rate varies with time as 

ip(t) oc eacp(— 7*) (4) 

from t = to the present galaxy age t — to- The parameter 7 is 
the inverse time-scale of star formation. For simplicity, all models 
are assumed for the moment to have fixed solar metallicity. 

Our main goal is to identify the model parameter that most 
simply accounts for the offset of normal, quiescent star-forming 
galaxies from the mean IRX-UV relation for starburst galaxies 
in Fig.0 To do this, we must explore the influence of each pri- 
mary model parameter (to, 7, tv and fi) on observable quanti- 
ties available for galaxies in our sample, i.e. the ultraviolet spec- 



tral slope, the ratio of far-infrared to ultraviolet luminosities and 
the 4000 A discontinuity. The influence of fy and /x on /3glx and 
L dustl 'J/fuv can be infe rred from the results presented in figure 2 
of ICharlot & FaM <200Cl) : these parameters, in particular the effec- 
tive dust absorption optical depth in the ambient ISM, fifv, control 
the position of a galaxy along the IRX-UV relation. The narrow- 
band spectral index D n (4000) is only weakly sensitive to variations 
in tv and /1 (Section 2). We find that the influence of the star forma- 
tion history on the IRX-UV relation through the parameters to and 
7 is best described in terms of the ratio of present to past-averaged 
star formation rates. This is defined by 

b = ^0) = to Mtp) 

W f ° dt'W)' 

For the expression of ip(t) in equation we have 

b _ 7ft) exp (-7*0) 
1 — cxp (—7ft)) 

For reference, a galaxy with constant star formation rate (7 = 0) 
has b = 1 at all ages. 2 

Fig-Ill summarizes me main results of our analysis. In Fig.|3ji, 
we show Ldu S t/-£/Fuv as a function of /3glx for sequences of 
models with different ages for four different time-scales of star for- 
mation, corresponding to 7 = 0.1, 0.5, 1.0 and 2.0 Gyr -1 , and for 
tv = 1.5 and fi = 0.3. The sequences are plotted as solid lines at 
ages when the ratio of present to past-averaged star formation rates 
is larger than b = 0.3 and as dashed lines at greater ages, down to 
b — 0.03. From b — 1 (age t — 0) to b = 0.3, the distance from 

2 Since the age to of the oldest stars in a galaxy is difficult to estimate 
observationally, it is often replaced b y some reference age, for exam ple the 
age of the Universe (e.g. Kennicutt 1983; Kennicutt, Tamblvn & Conadon 
1994). Such an approximation leads to larger b than obtained when using 
the definition given in equation j5J- 
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Figure 3. (a) Age sequences of models with star formation histories given 
by equation BT with 7 = 0.1, 0.5, 1.0 and 2.0Gyr _1 (as indicated) and 
for fy = 1.5 and fi = 0.3. The sequences are plotted as solid lines at 
ages when the ratio of present to past-averaged star formation rates is larger 
than b = 0.3 and as dashed lines at greater ages, down to b = 0.03. (b) 
Sequences of models with different dust content (i.e. different fy ) and /1 = 
0.3, for two ratios of present to past-averaged star formation rates, b = 
1.0 and 0.3 (as indicated), (c) Age sequences of models with 7 = 0.15 
(equation 0, which undergo an extra burst of star formation of duration 
10 s yr at the age of 6 Gyr involving 3, 10 and 30 per cent of the total mass 
of stars formed by the continuous model over 13.5 Gyr (as indicated). The 
dust parameters are fy = 1.5 and /i = 0.3. In all panels, the data and the 
mean starburst IRX-UV relation shown as a dotted line are the same as in 
Fig.Hl 



the mean starburst relation increases as b decreases along the model 
sequences, models with b w 0.3 occupying regions of the diagram 
populated by more quiescent galaxies. This behavior arises mainly 
from an increase in /3glx as b decreases. Models with b ~ 0.3 are 
dominated by stellar populations with intrinsically redder spectra 
and hence larger ultraviolet spectral slopes than starburst galaxies. 
We note that, for b > 0.3, all models in Fig.[3Ji follow a similar path 
in the IRX-UV diagram. For b < 0.3, the evolution in the IRX-UV 
diagram depends more strongly on the adopted time-scale of star 
formation. At fixed b, /3glx is larger in models with short time- 
scales of star formation (large 7), which contain a smaller propor- 
tion of ultraviolet-bright stars than models with extended star for- 
mation. In this b range, Z/dust/^ruv increases at late ages because 
of the contribution to dust heating by the optical radiation from old 
stars (Section 4). We note that models with small b 'returning' at 
late ages to the mean IRX-UV relation for starburst galaxies are 
considerably fainter than starburst galaxies and have intrinsically 
red optical spectra. 

The similarity for b > 0.3 between models with different 7 in 
the IRX-UV diagram prompts us to plot, in Fig.[3J>, I/dust /£fuv 
as a function of /3glx for sequences of models with different 



dust content (i.e. different fy), for two ratios of present to past- 
averaged star formation rates, b — 1.0 and 0.3. For simplicity, 
we assume that the fraction of the total effective absorption opti- 
cal depth contributed by the ambient ISM is fixed at the standard 
value fi — 0.3 in all models. The striking result of Fig.|3J) is that 
the model sequences with constant b are remarkably parallel to the 
mean IRX-UV relation for starburst galaxies, the b — 0.3 sequence 
being offset toward the region of the diagram populated by quies- 
cent star-forming galaxies. This is consistent with the observational 
evidence from Figs|2^ and b that, while dust content is the primary 
driver of the IRX-UV relation, the broadening of the relation for 
quiescent star-forming galaxies is most likely driven by star forma- 
tion history. This result does not depend on the specific combina- 
tion of 7 and to used to compute models with b = 0.3 in Fig.[5J). 
We find that, for b <C 0.3, sequences of models with different dust 
content are also parallel to the mean IRX-UV relation for starburst 
galaxies. However, as expected from Fig.|3^, the location of the se- 
quences in this b range depends sensitively on the choices of 7 and 
to (see also Section 4). 

Figs|3ji and b suggest that models with simple continuous 
star formation histories cannot easily reproduce the observations of 
galaxies lying at the largest perpendicular distances from the mean 
starburst IRX-UV relation in our sample, unless very short time- 
scales of star formation (7 ~ 2) and very small ratios of present 
to past-averaged star formation rates (b ~ 0.03) are invoked. Such 
parameters are not typical of quiescent star-forming galaxies. In 
Fig.[3j:, we show Ld us t / £fuv as a function of /3glx for models 
with a star formation time-scale parameter 7 = 0.15, more typical 
of a spiral galaxy, which undergo an extra burst of star formation of 
duration 10 8 yr at the age of 6 Gyr, for fy = 1.5 and fi — 0.3. We 
show the time evolution of these models for bursts representing 3, 
10 and 30 per cent of the total mass of stars formed by the contin- 
uous model over 13.5 Gyr. When the burst starts, each model first 
joins the mean IRX-UV relation for starburst galaxies. Shortly after 
the burst ends, Ldust and hence Ldust / £fuv drop rapidly because 
of the dispersal of the (optically thick) birth clouds of the last stars 
formed during the burst. Then, /9glx and I/dust /£fuv increase 
again as the ultraviolet light from these stars fades and reddens. 
A maximum in /?glx is reached about 0.5-1 Gyr after the burst 
ended. The subsequent evolution toward smaller /9glx, as the un- 
derlying continuous star formation again starts to dominate the pro- 
duction of the ultraviolet light, is remarkable for two reasons. First, 
it is very slow, implying that the models spend several gigayears 
in the region of the diagram populated by galaxies far away from 
the mean starburst IRX-UV relation. Second, the 6 parameter is 
only slightly lower than that of the undisturbed continuous model 
and hence typical of quiescent star-forming galaxies (note that this 
conserves the trend of decreasing b with increasing distance from 
the mean starburst relation). Thus, models including minor bursts 
of star formation on top of continuous star formation histories can 
account in a natural way for the observations of galaxies at large 
distances from the mean starburst relation in the IRX-UV diagram. 

The simple model described above provides a natural frame- 
work for understanding the physical origin of the IRX-UV diagram. 
As shown in Fig.|3| this model can account for the observational 
evidence from Section 2 that, while dust content is the primary 
driver of the IRX-UV relation for starburst galaxies, the broaden- 
ing of the relation for quiescent star-forming galaxies is most likely 
driven by star formation history, and in particular the 6 parameter. 
These results do not depend sensitively on the population synthesis 
code used for computing the production of starlight in the model. 
We have checked that adopting the PEGASE population synthesis 
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code iFio c & Rocca-Volm erangdl 19971 : version 2.0) instead of the 
iBruzual & Chariot! J2003T code has a negligible influence on the 
results of Fig.|3] 

It is worth noting that other effects, such as anisotropic emis- 
sion from galaxies and spatial distribution of the dust, can poten- 
tially affect the positions of galaxies in the IRX-UV diagram. The 
inclination angles available for 59 quiescent star-forming galaxies 
observed with /t/E+SDSS, ANS and OAO in our sample span the 
full range from to 90 degrees and do not exhibit any specific trend 
with perpendicular distance from the mean IRX-UV starburst rela- 
tion. Therefore, we do not expect orientation effects to play a major 
role in the b roadening of the IRX-UV relation for these galaxies. 
Charlo t"&FaTil 1 2000 ) show that, for some combinations of opti- 
cal properties and spatial distributions of the dust in their model, it 
is possible to produce I/dust/^FUv ratios lower than observed for 
starburst galaxies at large ultraviolet spectral slopes. However, the 
results in their figures 6 and 13 indicate that such models could not 
account for the low Ldust / £fuv ratios of quiescent star-forming 
galaxies with small /Sglx in our sample. 3 Changes in the spatial 
distribution of the dust also do not provide a natural explanation 
for the strong drop in Ha equivalent width with increasing distance 
from the mean starbu rst IRX-UV relation in Fig.[2f) (see figure 6c of 
ICharlot & F all 2000). We conclude that, although orientation and 
dust geometry may have minor effects on the positions of galaxies 
in the IRX-UV diagram, they do not seem responsible for the over- 
all broadening of the IRX-UV relation for quiescent star-forming 
galaxies in our sample. This broadening can be most simply under- 
stood as a sequence in the star formation histories of the galaxies. 



4 CONSTRAINTS ON THE ULTRAVIOLET 
ATTENUATION 

The results presented in the previous sections have important impli- 
cations for estimates of the ultraviolet attenuation Afuv and hence 
of star formation rates in galaxies. Meurer etalJ Jl999i) derive sim- 
ple formulae for estimating Afuv from either the ultraviolet spec- 
tral slope or the ratio of far-infrared to ultraviolet luminosities on 
the basis of the tight IRX-UV relation for starburst galaxies. The 
results of Sections 2 and 3 above suggest that these formulae can- 
not be straightforwardly applied to more quiescent galaxies. We 
can use the model described in Section 3 to estimate the accuracy 
to which Afuv can be constrained in galaxies of different types. A 
way to assess this is to investigate the relations between Afuv, 
/3glx and I/dust /£fuv for a library of models encompassing 
a wide range of physically plausible star formation histories and 
dust contents. We thus generate a library of Monte Carlo realiza- 
tions o f different star formation histories similar to that used by 
iKauffmann et alj <2003l) to interpret the D n (4000) and H5 A index 
strengths of a complete sample of over 10 J SDSS galaxies. 

The main requirement for the model library is to include broad 
enough ranges of star formation histories and dust contents to en- 
sure a wide coverage of the IRX-UV diagram. We have checked 
that the results presented below do not depend sensitively on the 
detailed parametrization of the models. For simplicity, we follow 



3 IWitt&Gordoi] 120001) also explore the influence of the optical properties 
and spatial distribution of the dust on the IRX-UV relation using a code 
of radiative transfer with spherical symmetry. None of the models they in- 
vestigate produces idust/^FUV ratios lower than observed for starburst 
galaxies at fixed ultraviolet spectral slope (see their figure 12, mislabeled 
figure 13 in the print version). 



IKauffmann et alj l2003li and parametrize each star formation his- 
tory in terms of two components: an underlying continuous model 
with a star formation law given by equation 0, and random bursts 
superimposed on this continuous model. 4 For the purpose of the 
present analysis, we retain only models with ratios of present to 
past-averaged star formation rates larger than 0.03. For the attenua- 
tion by dust, we adopt a broad distribution of fv with a peak around 
1.5 and a broad distribution in /j, with a peak around 0.3. This is 
consistent with the dust properties derived from the Ha/H/3 ratios 
and broad-band optical magnitud es of the above sample of 10 5 
SDSS galaxies iBrinchmann et alJl2003l: IKauffmann et alj|2003h . 
We distribute our models logarithmically in metallicity from 0.1 to 
2 times solar. Our final library consists of 95,000 different models. 

Fig-El shows ^dust/i/Fuv as a function of /3glx for these 
models. In the upper panel, we have divided the IRX-UV plane 
into bins and coded each bin with a different colour to reflect the 
average ratio of present to past-averaged star formation rates of the 
simulated galaxies that fall into it. As found in Section 3 above, se- 
quences of galaxies with different dust contents corresponding to 
different b lie parallel to the mean IRX-UV relation for starburst 
galaxies. In the lower panel of Fig.|4| we have coded each bin in 
the IRX-UV plane with a different colour to reflect the average ul- 
traviolet attenuation Afuv for the galaxies that fall into it. It is 
remarkable that sequences of models with different star formation 
histories corresponding to different Afuv are nearly perpendicular 
to the y-axis in this diagram. This implies that Ldust/^FUV is a 
much better estimator of Afuv than /3glx in galaxies with differ- 
ent star formation histories. 

To quantify the accuracy to which Afuv can be estimated 
from either Ldust /-^fuv or /3glx, we have explored functional 
fits of Afuv in terms of these parameters using the models in our 
library. We find that the dependence of Afuv on Ld us t / Lfuv is 
well approximated by two linear relations with a change in slope at 
some characteristic turn-around point (corresponding to Afuv ~ 
1). This can be represented by an expression of the form 



Ar 



(Ao + Ai)(x- xo) 



+ 



2 

Ao)cr 



+ 



In 



cosh 



(Hr) 



+ A 2 



(7) 



where x = log(Ld us t/£Fuv) and Ao (the slope of the relation at 
small x), Ai (the slope of the relation at large x), Ai (the normal- 
ization constant), xo (the turn-around point) and ao (the smooth- 
ness of turn-around) are fitted parameters. The accuracy of the 
Afuv estimates obtained from this formula depends on whether 
or not independent information is available on the ratio of present 
to past-averaged star formation rates. If no information is avail- 
able on b, equation (7J with A = 0.39, Ai = 1.86, A 2 = 0.75, 
xo = 0.26 and ao = 0.42 provides estimates of Afuv with a typ- 
ical (i.e. rms) relative uncertainty of 20 per cent for all the models 
in the library. 

4 Specifically, we take the galaxy age to to be distributed uniformly over 
the interval from 1.5 to 13.5 Gyr and the star formation time-scale parameter 
7 over the interval from to 1 Gyr~ 1 . Random bursts occur with equal 
probability at all times until to ■ They are parametrized in terms of the ratio 
between the mass of stars formed in the burst and the total mass of stars 
formed by the continuous model over the time to ■ The ratio is taken to be 
distributed logarithmically from 0.03 to 4.0. During a burst, stars form at a 
constant rate for a time distributed uniformly in the range 3 X 10 7 -3 X 10 8 
years. The burst probability is set so that 50 per cent of the galaxies in the 
librar y have experienced a burst in the past 2 Gyr (see Kauffmann et al. 
2003). 
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Table 1. Parameters of the fits of equation Q f° r different ranges of present to past-averaged star formation 
rates b. Also indicated for reference are the corresponding ranges in specific star formation rate i/> (to ) / M t 
(a quantity more loosely connected to the star formation history but easier to constrain observationally). The 
last row shows the results when no restriction is applied on b and hence ip(to)/ M* . Also listed is the typical 
(i.e. rms) relative uncertainty s rc i in the estimated j4_fuv for the library of models described in Section 4. 



Range in log b 


Range in log { [^{to)/M»] (vr" 1 )} 


A 


Ai 


A 2 


x 




•^rcl 


[-1.5, -1.0] 


[-11.3, -10.6] 


0.27 


2.02 


1.01 


0.69 


0.91 


0.15 


[-1.0, -0.5] 


[-10.6, -10.0] 


0.34 


2.07 


0.90 


0.44 


0.62 


0.10 


[-0.5, +0.0] 


[-10.0, -9.3] 


0.35 


2.24 


0.92 


0.36 


0.60 


0.06 


[+0.0, +2.0] 


[ -9.3, -7.5] 


0.39 


2.37 


0.96 


0.34 


0.59 


0.09 


Full range 


Full range 


0.39 


1.86 


0.75 


0.26 


0.42 


0.19 




Figure 4. Ratio of far-infrared to ultraviolet luminosities plotted against 
ultraviolet spectral slope (the TRX-UV diagram') for the model library de- 
scribed in Section 4. Upper panel: the IRX-UV plane has been binned and 
coded to reflect the average ratio of present to past-averaged star forma- 
tion rates of the galaxies that fall into each bin. Lower panel: the IRX-UV 
plane has been binned and coded to reflect the average ultraviolet attenua- 
tion Aprjv of the galaxies that fall into each bin. In each panel, the solid 
line shows the mean starburst IRX-UV relation from Fig. HI 

More accurate estimates of Afuv can be obtained when even 
rough information is available on b. This is illustrated in Tabled 
in which we list the parameters of the fits and the accuracy of the 
Afuv estimates obtained from equation Q f° r different ranges in 
log b. For reference, we also indicate the corresponding ranges in 
'specific star formation rate' ip(to) / M* (i.e. the current star forma- 
tion rate per unit stellar mass), a quantity more loosely connected 
to the star formation history but easier to constrain observationally. 
The typical relative uncertainty in the Afuv estimates decreases 



from 15 per cent for the most quiescent galaxies (log b < —1) to 
only 7 per cent for starburst galaxies (log b > 0). The reason for 
this is that, in starburst galaxies, the heating of the dust and the 
production of the ultraviolet radiation are dominated by the same 
stars. Thus idust/iFUV is intimately related to Afuv- In quies- 
cent galaxies, Ld us t /I/fuv is affected by dust heating by the op- 
tical radiation from old stars with little ultraviolet emission, which 
depends on the history of star formation at fixed Afuv. For ref- 
erence, the contribution to dust heating by stars older than 0.3 Gyr 
ranges from typically less than 10 per cent for b > 1 to about 80 
per cent for b ~ 0.03 in the models of Fig.|4| This is the reason 
why, when using the fitting formulae derived from TableQ the same 
Afuv is obtained for l arger Ldust /Lfuv when b is lower. A sim- 
ilar result was noted bv lBuat et al.Hl99H) in their analysis of a far- 
infrared selected sample of nearby star-forming galaxies. We con- 
clude that, if rough constraints are available on the ratio of present 
to past-averaged star formation rates, reliable estimates of Afuv 
can be derived from idust/ipuv- We note that, for the models in 
our library, log b can be roughly estimated from D n (4000) as 

log b = -2.52 + 40.81 exp [-2.24 D n (4000)] , (8) 

with a dispersion of about 0.3 at fixed D n (4000). The corre- 
sponding dependence of the specific star formation rate V>(tn) /M, 
on D n (4000) is in ex cellent agreement with that found by 
iBrinchmann et alJ <2003l) for the sample of 10 5 SDSS galaxies 
mentioned above. 

In some cases, the far-infrared luminosity may not be avail- 
able to estimate Afuv (for example, for high-redshift galaxies for 
which only the rest-frame ultraviolet emission can be observed). In 
such cases, it is important to know the accuracy to which Afuv 
can be estimated based solely on the ultraviolet spectral slope. The 
wide range in Afuv spanned by models with fixed /3glx in Fig.|4] 
indicates that /3glx is not a precise estimator of Afuv in galaxies 
with diff erent star fo rmation histories. This is consistent with the 
results o f lBeljfcOOA . In fact, we did not succeed in finding a func- 
tion of /?glx that provided reasonably accurate estimates of Afuv 
for the models in our library. Significant improvement is obtained 
when including the ratio of present to past-averaged star formation 
rates as a variable of the fit. In particular, we find that the depen- 
dence of Afuv on /3glx and b is reasonably well approximated 
by 

Afuv = 3.87 + 1.87 ( /3 G lx + 0.40 log b) . (9) 

This formula provides estimates of Afuv for which the accuracy 
is an increasing function of b. For b > 0.3 (corresponding roughly 
to types later than Sb), the typical (i.e. rms) absolute uncertainty 
in the Afuv estimates is only 0.32. For lower b, the uncertainty 
is about one magnitude. The reason for the large uncertainty in 
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j4fuv at low b is that, in galaxies with little current star formation, 
the ultraviolet spectral slope depends strongly on both dust and the 
history of star formation. We note that, for b » 5, the expression 
of Afuv given in equation {5J is similar to the non-6-dependent 
expression for A-pw as a function of ultraviolet spectral slope de- 
rived by Meurer et al. (1999, their equation 1 1) from the analysis of 
an ultraviolet-selected sample of nearby starburst galaxies. In the 
absence of far-infrared observations, therefore, j4fuv can still be 
reasonably well estimated from the ultraviolet spectral slope alone, 
provided that some constraints are available on the ratio of present 
to past-averaged star formation rates. 



5 SUMMARY AND CONCLUSION 

We have compiled a sample of 115 nearby, non-Seyfert galaxies 
drawn from ultraviolet surveys spanning a wide range of star for- 
mation activities, from starburst to nearly dormant. Based on this 
sample, we showed that, at fixed ultraviolet spectral slope ( /3glx), 
normal star-forming galaxies have systematically lower ratio of to- 
tal far-infrared to ultraviolet luminosities (Ldust/^FUv) than star- 
burst galaxies. The tight ness of the cor relation between /3glx and 
Ldvist/ L-pvw, shown by Meurer et"al] \ 19991) to provide valuable 
constraints on the ultraviolet attenuation in starburst galaxies, does 
not extend therefore to normal, more quiescent star-forming galax- 
ies. We showed that, for the galaxies in our sample, the strengths 
of the age-sensitive 4000 A spectral discontinuity and Ha emission 
equivalent width correlate well with distance from the mean rela- 
tion for starburst galaxies in the IRX-UV diagram, while there is 
little or no correlation between the dust-sensitive Ha/H/3 ratio and 
this distance. This is strong observational evidence that, while dust 
content is the primary driver of the IRX-UV relation for starburst 
galaxies, the broadening of the relation for quiescent star-forming 
galaxies is most likely driven by star formation history. 

We hav e used the simple bu t physically realistic attenua- 
tion model of Ich arlot & Fall 1 2000), in combination with the new 
iBruzual & Chariot! 120031) population synthesis code, to investigate 
the physical origin of the IRX-UV relation for galaxies with differ- 
ent star formation histories. We confirm that, for starburst galaxies, 
the tight IRX-UV relation can be understood most simply as a se- 
quence in overall dust content. In addition, within the framework 
of this model, the second dimension of the IRX-UV relation can 
be understood most simply as a sequence in the ratio of present to 
past-averaged star formation rates. 

The dependence of the IRX-UV relation on both star forma- 
tion history and dust content implies that the ultraviolet attenu- 
ation Afuv cannot be derived from I/dust /£fuv and /3glx as 
straightforwardly as expected from the analysis of starburst galax- 
ies alone. We have used a library of Monte Carlo realizations of 
galaxies with different star formation histories and dust contents to 
quantify the accuracy to which Afuv can be estimated from ei- 
ther I/dust /Lfvv or /3glx- We find that L dust /L FV v is a better 
estimator of Afjjv than /3glx in galaxies with different star for- 
mation histories. We provide simple formulae for estimating Afuv 
as a function of either Ldust/^FUV or /3glx and show that the 
accuracy of these estimates can be significantly improved if some 
constraints are available on the ratio of present to past-averaged 
star formation rates (for example, based on the observed 4000 A 
discontinuity). The IRX-UV diagram, therefore, is a valuable diag- 
nostic of the attenuation by dust and hence the star formation rates 
of galaxies in a wide range of star formation activities. This is espe- 
cially important for applications to studies of high-redshift galax- 
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ies, for which the rest-frame optical and infrared emission cannot 
be observed at present. 
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Table Al. Effective wavelengths and FWHMs of the ultraviolet filters used 
on-board different satellites. 



Satellite 


A e (A) 


FWHM (A) 


GALEX (FUV) 


1520 


264 


GALEX (NUV) 


2310 


795 


ANS (1) 


1550 


149 


ANS (2) 


1950 


149 


ANS (3) 


2200 


200 


ANS (4) 


2500 


150 


OAO (1) 


1550 


240 


OAO (2) 


1910 


260 


OAO (3) 


2460 


380 


UIT(Bl) 


1521 


354 


UIT (Al) 


2488 


1147 



APPENDIX A: CALIBRATION OF ULTRAVIOLET 
SPECTRAL SLOPES AND FLUXES 

Here, we use the sample of 110 high-signal-to-noise IUE galaxies 
with both far- and near-ultraviolet spectra compiled in Section 2. 1 
to calibrate the conversions between /3glx (equation^ and the 
ultraviolet spectral slopes estimated from observations with differ- 
ent satellites. The minimum signal-to-noise ratio per pixel 5 of 3 
required for this sample ensures that /3glx can be measured to 
an accuracy of better than 0.3. For reference, the median error in 

/3glx for the 1 10 galaxies in our sample is around 0.08. 

We compute the ultraviolet spectral slopes using vari- 
ous definitions and compare the results to /3glx- We first 
consider the ultraviolet spectral slope fro m the definition of 
ICalzetti. Kinnev & Storchi-BergmannHl994l) . which we denote by 

/3cks- This is obtained from a power-law fit of the form fx oc X 13 
(where fx is the flux per unit wavelength) to the spectrum in ten 
continuum bands that avoid strong stellar and interstellar absorp- 
tion features (in particular, the 2175 A dust feature) in the range 
1268 ^ A 2580 A. We also consider the ultraviolet slopes de- 
rived from photometric observations with ANS, OAO and UIT. We 
compute the corresponding Pans, Poao and Pun from contin- 
uum fits of the form fx oc A^ to the broadband spectra obtained 
by convolving the IUE spectra with the filter response functions of 
these satellites (Table lAfl , 

Fig.|AJJshows /3cks, /3ANS,and (3oao and /Juit as a func- 
tion of /3glx for this sample. It can be seen that Pens, Pans 
and Poao correlate extremely well with /3glx- The correlation is 
worse for /3uit, mainly because of the wider wavelength range of 
the near-ultraviolet filter of UIT, which extends redward to about 
3060 A (Table|A3. The mean linear relations between the ultravi- 
olet slopes obtained from observations with different satellites and 
the reference ultraviolet slope /3glx are drawn as solid lines in 
Fig. IA II and are indicated at the bottom of each panel along with 
the lcr standard deviation. The tight correlations in Figs. l(a)-l(c) 
allow us to estimate /9glx reliably from Pcks Pans and Poao • 
Since the relation is significantly worse for /Suit, we do not use 
UIT data in the analysis described in the main text of the paper. 



The ultraviolet spectral slope /3glx can also be estimated re- 
liably for galaxies with only far-ultraviolet IUE spectra. The reason 
for this is that the wavelength range of the far-ultraviolet spectro- 
graph on-board IUE encompasses 9 out of the 10 spectral windows 
originally used by Calzetti et alj J 19941) to define /3cks- In fact, 
for the 110 galaxies in the sample considered here, we find that 
the ultraviolet slope P'oks obtained by fitting a power law to the 
far-ultraviolet spectra in these 9 continuum bands correlates with 
/3glx as tightly as do Pans and Poao ■ The mean linear relation 
is 

Pglx = ( P'cks ~ 0.38)/1.05 , (Al) 

with a standard deviation a — 0.23. 

We also adopt a common definition of the ultraviolet flux for 
galaxies observed with different satellites, Jfuv = Afuv/fuv- 
For those galaxies with IUE spectra, we compute the mean ultravi- 
olet flux density / FUV using the response function of the GALEX 
far-ultraviolet filter (Afuv = 1520 A). For those galaxies ob- 
served with OAO and ANS, we compute / FUV by extrapolating the 
broadband spectrum from 1550 A (corresponding to the effective 
wavelength of the far-ultraviolet filters on-board both satellites) to 
1520 A, using the ultraviolet slope ( Poao or Pans) measured for 
each galaxy. Given the narrow wavelength range over which the 
extrapolation is performed, this procedure does not introduce any 
significant uncertainties in / FUV and hence Ffuv- 



5 This is the average signal-to-noise ratio per pixel in 100 A-wide contin- 
uum windows free of e mission features cen tred at rest wavelengths 1450, 
1700, 2300 and 2700 A iKinnev et alll99lh . 
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Figure Al. Ultraviolet spectral slopes computed in the photometric bands observed by different satellites plotted as a function of /3glx. for the 110 IUE 
galaxies with both far- and near-ultraviolet spectra in the sample of Section 2.1. In each panel, the inset diagram shows the distribution of the difference 
between the true value of /?glx and that obtained from the mean linear relation (indicated at the bottom and drawn as a solid line) for this sample, along with 
the lcr standard deviation. 



